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INTRODUCTION 


It has been reported that annular pressure seals in centrifugal pumps 

and Vr a ? lnf i U£nCe on the lateral vibrations of the shaft systems (ref^ 1 
and 2 ). In order to analyze the dynamic properties of the fluid films in the 
annular seals, it is necessary to know the boundary conditions of pressure 
such as pressure at the both ends of the seal. It is also important to evalu- 

! the res ytance coefficient of the flow in the seal both in the 

axial and the circumferential directions. Unfortunately there are few data 

s^al is :omp°iiLle U d a ' e ^ ValU6S ’ when the configuration of the 


In this study two series of experiments, 
are carried out. In the static tests, pressure 
seals are measured under several conditions of 
pressure difference and the speed of the rotor 
the exit loss coefficients and the resistance 
the dynamic analysis. 


static tests and dynamic tests, 
distributions in the eccentric 
the eccentricity, the axial 
to investigate the inlet and 
coefficients, which will help 


Then the dynamic forces at the 
measured. The dynamic properties of 
lated from the test results. 


seals caused by a vibrating rotor are 
the fluid films in the seals are calcu- 
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NOMENCLATURE 
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The following nomenclature is used in the paper: 


damping coefficient, N-s/m 
seal diameter, m 
eccentricity ratio ( =6/G ) 

static force on an eccentric rotor, N ^ 

non-dimensional force on the rotor ( =f / (0 . 5p • v m *d- 1) ) 

dynamic force on the stator, N 

seal clearance when concentric, m 

stiffness coefficient, N/m 

seal width, m 

inertia coefficient, N»s /m 
rotational speed, rps 
pressure, N/m 

non-dimensional pressure ( =(p-p^ ) / (p, a - 

non-dimensional maximum pressure difference in a section perpendicu ar 

to the axis ( = (P^-P^x-x / ) 

leakage flow rate, m 3 /s 

axial Reynolds number ( =v m *G/v ) 

rotational Reynolds number ( =n- d- N * G/v ) 

axial velocity component of the flow in the seal, m/s 

mean axial velocity ( =Q/(TT-d*G) ) , m/s 

co-ordinate in the axial direction ( x=0 at the inlet of the seal ) , m 
rotor displacement in the vertical direction, m 
rotor displacement in the horizontal direction, m 


eccentricity, m 

co-ordinate in the circumferential direction ( 


0=0 at the top of the 


seal ) 2 

resistance coefficient ( dp=X* (dx/2G) - (p ,, v rl /2) ) -oaf. z ^ 

resistance coefficient proposed by Yamada ( =0 . 26R a jl+( 7 /8) • (Rr /2R a )) , 

kinematic viscosity, m 2 /s 2 . 

inlet loss coefficient ( = (p, -p 2 ) / (0 . 5p * v m ) -X* x 2 / ( 2G) -1 ) 
exit loss coefficient ( =l-(p 6 - P 4 ) / (0 . 5p* v£) -X’ ( l-x 4 ) / ( G) ) 


density, kg/m 

vibration angular velocity , rad/ s 


Subscripts : 

0 amplitude at the frequency given by the vibrator 
1-6 value at the measuring position 1-6 

p _3 value corresponding to the vibration angular velocity co, Wj 
X component in the vertical direction 

Y component in the horizontal direction 

XX XY, YX, YY used in the dynamic properties of the fluid film; The irs 
subscript denotes the direction of the force. The second subscript 
denotes the direction of the displacement, velocity or acceleration. 
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STATIC TEST APPARATUS 


The static test apparatus is shown in figure 1. Static pressure was 
measured at the pressure holes distributed at six axial positions and eight 
circumferential positions with pressure transducers of strain gauge type. 

The eccentricity of the rotor was measured with displacement sensors of eddy 
current type at two axial positions and two circumferential positions 90 
degree apart from each other. The radial load on the bearings in the both 
vertical and horizontal directions were also measured with load transducers 
of strain gauge type. The static pressure in the high and low pressure cham- 
bers were measured as well as the flow rate through the seals. 

Tests were performed under conditions of the axial Reynolds number R a 
between 1100 and 3500 and the rotational Reynolds number R r between 0 and 
4800. 


TEST SEALS 


The configurations of test seals are shown in figure 2. The seals have 
relatively small ratio of width to diameter. Figure 2 also shows the axial 
positions of the pressure holes 1 to 6 and the axial positions of the dis- 
placement sensors A and B. 


TEST RESULTS 


Figure 3 shows a typical pressure distribution measured in the test seal 
I . It is noted that the maximum pressue in a section perpendicular to the 
axis was measured near the circumferential position where the clearance was 
minimum throughout the whole width of the seal. It was the same regardless 
of the pressure difference or the rotor speed, provided that the ratio of the 
width to diameter of _the seal was small and that the seal had grooves iust 
like the seals I and Ui. shown in figure 2. 


It should also be noted that even at the position 5 behind the seal exit 
considerable pressure difference in the circumferential direction still re- 
mained. This result implies that it is not appropriate to the analysis of the 
flow to assume that the pressure at the exit of the seal is uniform and equal 
to the pressure in the low pressure chamber. 


Figure 4 shows the maximum pressure difference "AoV in the section perpen- 
dicular to the axis. In the case of the flat seal ( test seal E ) 7Tp de- 
creased linearly in the axial direction. On the other hand, in the cTse of 

the grooved seals ( test seals I and I ) 7^ had the minimum value and then 
increased again. 
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Figure 5 shows the influence of the inclination between the axes of the 
rotor and the stator on the pressure distributions in the seal. The effect of 
the inclination can not be neglected, especially when the eccentricity is 
large. 

Figures 6 and 7 show the inlet and the exit loss coefficients respective- 
ly calculated from the pressure distributions. The inlet loss coefficient 
was much smaller than the conventional value, 0.5. The exit loss coefficient 
Cout was also considerably smaller than the conventional value, 1.0. Since the 
pressure at the exit is expressed by the equation 

p(0)=p 6 -(l-^.p-v(0) 1 /2 (1) 

small gives a reasonable explanation to the pressure distribution there. 
Acording to the experimental results by Stampa ( ref. 3 ), inlet and exit loss 
coefficients are strongly affected by the axial and the rotational Reynolds 
numbers. Within the range of the Reynolds numbers of the experiments presented 
here, the rotational Reynolds number did not have much influence on the inlet 
loss' coefficient, though it had some influence on the exit loss coefficient. 

Figure 8 shows the resistance coefficient A calculated from the pressure 
distribution. It was normalized by A 0 , which was calculated by the equation 
proposed by Yamada for concentric flat seals ( ref. 4 ). The experimental 
resistance coefficient A for the flat seal agreed well with the coefficient 
calculated by Yamada ’s equation with consideration for the effect of eccen- 
tricity analyzed by Tao and Donovan ( ref. 5 ) . 

9 shows the non-dimensional radial forces calculated from the 
measured pressure distributions. In the case of the test seals I and BL , the 
radial forces did not depend much on the Reynolds numbers in these experiments. 

On the other hand in the case of the test seal E , the ratio of the rotational 

Reynolds number to the axial one had considerable influence on the forces. It 

is noted that the flow resistance in the circumferential direction is large 

in the flat seal and the rotor revolution has considerable effect on the flow 
in the eccentric flat seals. 


DYNAMIC TEST APPARATUS 


Figure 10 shows the dynamic test apparatus. The shaft was vibrated verti- 
cally by a hydraulic vibrator. Vertical and horizontal relative vibrations 
between the rotor and the stator of the test seal were measured at two axial 
positions with displacement sensors of eddy current type. Vertical and hori- 
zontal forces acting on the stator were measured at two axial positions with 
load sensors of strain gauge type. 


TEST DATA PROCESSING 

Test data were recorded on magnetic tapes and analyzed with a computer. 
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First the analogue signals were converted into digital signals with a high 
speed AD converter. Then the digital signals were processed by the digital 
filter method and the amplitudes and the phases at the frequency given by the 
vibrator were calculated. Vibration amplitudes at two axial positions were 
averaged to be used in the farther processing. Load amplitudes at two axial 
positions were added to each other to obtain the total values. The acceleration 
of the casing was also measured and was taken into account to correct the load 
amplitudes. The amplitudes and the phases were transformed into the real and 
the imaginary parts of complex expressions. 


The relations between the dynamic properties of the fluid film and the 
hydraulic forces acting on the stator are expressed by the following equations 

m xx^ + m xy Y " + c AX X + c xf Y + k xx X + k x y Y = F x 

( 2 ) 

n \x X + m Yr Y + c^X + c yY Y + k^X + k yr Y = F Y 

The vibrations of the rotor and the forces on the stator may be expressed in 

the complex number as follows. 

X = X 0 e'^ = { R e (X 0 ) + il w (x 0 )} e |lA,t 

Y = Y c e' wC = -jR e (Y t ,) + iWYc)} e ,u>L 

F X = F Xo e iwt = jR e .(F X0 ) + iI*(F xc )} e"^ ^ 

F r = = jR e (F Y0 ) + iI m (F Yo )} e'^ 


Since only four equations are obtained for twelve unknowns m^ , m XY , m YX , m YY , 
c xx > , c Yx , c Yy -, k xx , k xr , k YX and k yy for each vibration frequency, three 

experiments at different vibration frequencies are necessary to obtain full 
equations. The final equations are expressed in the matrix form as follows. 


ReCXci) R* (Yi, ) -on, I m (Xo.) 
Re(XoJ Kc>a.) -ahLnCx*) 

Re (Xm) Re (^cy) -lib I m (Xcj) 

I/nCXc,) Ln (Yc , ) W t Re(Xo,) 

f**t 1« (Y^) ^2 Re (^ 02 ) 

Im(Y CJ ) oyReCXw) 


W| I M (Y c , ) 
-UJilw (Y w ) 
-00,1*^) 
Uj Ra (Yci ) 

Re (Y^) 
Re (Y^) 


-u?Re(X cl ) -u \ 2 R e (Yc,) 
-wi-Re (X«) -c^Re(Y a ) 

-U^RdCXrt) -^ReCYo,) 
-W^LnCXo.) 

I* ( Xj X ) —(ii* Im (Y ci ) 
«(\i) -<I«(Y W ) 


k xx 


Re(F«,, ) 

kxy 


Re (F^J 



Re(F^) 

C *Y 


I„(F,w,) 

m *x 


I*(F>,z, z ) 

m xr 


I m (FxOj) 


(4) 


Only six equations for the equilibrium of the forces in the X direction are 
shown here. The equations in the Y direction are obtained by replacing the 
unknown vector and the vector in the right hand side. 


TEST SEAL AND TEST RESULTS 

A test seal similar to the test seal I in figure 2 was set in the appa- 
ratus and tested. 

Figure 11 shows typical test data. Six vibration frequencies from 10 to 
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40 Hz were selected. The vibration amplitude R*(X C ) was set around 12 pm ( 1 L 
of the seal clearance ) and I m (X 0 ) and Y a were very small regardless of the 
vibration frequency . 

Re (Fxo ) /Rfc(Xe) decreased parabolically with the increase in the^ vibra- 
tion frequency because in these experiments Re.(F*c ) /Re(X c ) ^ “ ^ m xx • 

On the other hand I m (F* c ) /Re (Xc) increased linearly with the vibration fre- 
quency because I*(F X0 ) /R e(X 0 ) = (jJ-c xa . 

Figure 12 shows some results of dynamic properties of the fluid film in 
the seal. The pressure difference across the seal was maintained nearly con- 
stant. The time-averaged eccentricity of the rotor changed with the rotational 
speed. Damping coefficient c xx and inertia coefficient m xX did not change much 
with the speed or the eccentricity within the experimental range of the param- 
eters. Stiffness coefficient k xx changed considerably with the speed. From the 
results of the static test ( fig. 9 ) this change in k xx is considered to be 
the effect of the eccentricity rather than the effect of the speed. 


CONCLUDING REMARKS 


An experimental study was performed with static and dynamic test appa- 
ratus to investigate the hydraulic forces caused by annular pressure seals . 

The following conclusions are deduced. 

1. The pressure at the exit of the seal had considerable variation in the 
circumferential direction. 

2. The inclination between the axes of the rotor and the stator had much 
influence on the pressure distribution in the seal. 

3. The inlet loss coefficient was much smaller than the conventional value, 
0.5. It was affected by the axial Reynolds number. 

4. The exit loss coefficient was smaller than the conventional value, 1.0. 
This result gives a reasonable explanation to the conclusion 1. 

5. The experimental resistance coefficient for the flat seal agreed well 
with the coefficient calculated by Yamada's equation with consideration for 
the effect of eccentricity analyzed by Tao and Donovan. 

6. The hydraulic forces caused by the eccentricity of the rotor had strong 
non-linearity to the eccentricity. 

7. The results of the dynamic test showed that the damping coefficient Cx* 
and the inertia coefficient m xx were not affected much by the rotational speed 
or the eccentricity of the rotor, though the stiffness coefficient k xX seemed 
to be influenced by the eccentricity. 
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PRESSURE PORTS 



TEST SEAL I 


TEST SEAL II 


TEST SEAL III 
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Figure 2. Test seals 
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Figure 3. Pressure distribution in an eccentric seal 



Figure 4 . Change in circumferential pressure difference 

in the axial direction 
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INLET LOSS 
COEFFICIENT £in 



TEST SEAL I 

N = 16.7 rps 
P| ~ P6 =0.27 MP a 

a : inclination 


Figure 5. Pressure distribution in the axial direction 

-Influence of inclination between rotor and stator axes 




Figure 6. Inlet loss coefficient 


Figure 7. Exit loss coefficient 
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Figure 10. Experimental apparatus for annular seal dynamic properties 
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Figure 1 1 . An example of experimental data 






